Light Filaments Without Self Guiding 
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An examination of the propagation of intense 200 fs pulses in water reveals light filaments not 
sustained by the balance between Kerr-induced self-focusing and plasma-induced defocusing. Their 
appearance is interpreted as the consequence of a spontaneous reshaping of the wave packet form a 
gaussian into a conical wave, driven by the requirement of maximum localization, minimum losses 
and stationarity in the presence of non-linear absorption. 
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Since its discovery by Braun et al. Q in 1995, the 
spontaneous filament formation accompanying intense 
fs-pulse propagation in air has received rapidly increas- 
ing attention, both for the generation of coherent, soft 
X-rays |2|, IR |3]> as well as sub-terahertz m radiation, 
and for remote sensing in the atmosphere 5]. More re- 
cently research has been extended to the case of filament 
formation in condensed matter: namely, in fused silica Q 
and water 0, , mainly in connection with fundamental 
investigation in the soliton field |9|, . In spite of the dif- 
ferent power and length scales that the process exhibits, 
two key features emerge, which are substantially the same 
in all media investigated: (i) for powers well exceeding 
the critical value for continuous- wave (CW) self-focusing, 
a light filament appears and propagates in the absence 
of diffraction or optical breakdown for several to many 
diffraction lengths; and (ii) the filament only contains a 
small fraction of the total beam power, limited by the so- 
called "intensity clamping effect" , while the residual 
excess remains in the filament periphery with no appar- 
ent trapping. 

Two major approaches have been proposed for the 
modeling, description and understanding of the under- 
lying physics: the first interprets the filament as a 
genuine soliton-like, self-guided beam, whose station- 
arity is supported by the dynamical balance between 
Kerr-induced self-focusing and plasma-induced defocus- 
ing, where plasma formation is due to multi-photon ab- 
sorption (MPA) . In this scenario the role of the non- 
trapped radiation is marginal. In fact, although compen- 
sation of MPA-induced power losses by external radiation 
has been addressed experimentally pjj, this refilling ef- 
fect was only interpreted as perturbative, i.e., useful for 
prolonging the filament life while not being structural 
to its existence. The second approach, in contrast, only 
considers the filament as an optical illusion related to the 
use of time-integrated detection. The models, based on 
moving focus 14] or dynamic spatial replenishment [l5| . 
interpret the filament as continuously absorbed and re- 
generated by subsequent focusing of different temporal 
slices of the pulse. 
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FIG. 1: Profiles of (a) clipped, (b) stopped and (c) free fil- 
aments for 20 mm< z < 40 mm. The input peak power is 
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In this letter we demonstrate how, for the case inves- 
tigated of 200 fs filaments in water, energy refilling by 
the surrounding radiation is so important that one can 
by no means consider the filament in terms of a genuine, 
self-guided beam. Our results indicate that filament re- 
generation can well be described within the framework of 
a purely spatial process, with the fake self-guided beam 
as resulting from subsequent focusing of different spa- 
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tial coronas of the beam. The key picture we propose 
as a physical interpretation of the entire dynamics ad- 
dresses the spontaneous transformation of a Gaussian 
into a Bessel-type beam, driven by the requirements of 
minimum (non-linear) losses, maximum stationarity and 
maximum localization. 

The criterion we adopted for distinguishing between a 
true (e.g., self-guided) filament and a fake (e.g., strongly 
refilled) is based on the simple assumption that a gen- 
uine filament should survive (for at least 2-3 diffraction 
lengths) after transmission through a pinhole that re- 
moves the non-trapped part of the radiation and trans- 
mits most of the filament energy content. The Kerr 
medium that we have chosen is water; with respect to 
solid-state media, it has the key advantages of being free 
from damage problems, of allowing easy and continuous 
scanning of sample length, and of permitting the inser- 
tion of pinholes along the beam path inside the non-linear 
media; in comparison with air, water has much less tur- 
bulence, which results in higher beam-pointing stability 
and so allows the use of pinholes of comparable size to 
the filament. We note that this was not the case for the 
experiment in air quoted above 0, where a very large 
(e.g., 1mm) pinhole was used. 

The experiment was performed by launching a 527 nm, 
~3|uJ, 200 fs, spatially filtered beam with an approxi- 
mately 0.1mm FWHM waist located at the input facet 
of a water-filled cuvette, and then monitoring the output- 
beam fluence profile. The wave packet launched was 
provided by an SHG-compressed, CPA Nd:glass laser 
(TWINKLE, Light Conversion Ltd.), operated at a 33 Hz 
repetition rate. The cuvette was made of 1 mm thick, 
syringe-shaped quartz, which allowed continuous tuning 
of the sample length in the range z = 5-40 mm. The out- 
put beam was imaged onto a CCD camera (8-bit Pulnix 
TM-6CN) by an / = +50 mm achromatic objective, with 
8x magnification. The results highlighted the appear- 
ance of a single filament with almost constant ~ 20 jam 
FWHM diameter in the z = 15-40 mm range investi- 
gated. 

Figure GJa) shows the effect of inserting a 55 |im- 
diameter pinhole in the water cuvette at z = 20 mm. 
The energy transmitted by the pinhole was 20% of the 
total incident energy. Although the transmitted beam 
attempted to focus after a short distance, (e.g., within 
one diffraction length; see the plot at z = 22 mm), the 
filament survived no further and we observed a rapid de- 
cay with divergence ~ 2 x larger than that of a Gaussian 
beam of the same FWHM diameter. It is necessary to 
double check whether, in the absence of pinhole, the fil- 
ament survives just because it is strongly supported by 
energy refilling from the outside beam. To do this, we 
also performed the complementary experiment, by block- 
ing the central spike with the sole aid of a 55 \xm beam 
stopper, printed on a 100 urn-thick BK7 glass plate and 
inserted into the beam path at z = 20 mm. As depicted in 



Fig- Elk), a central spike of the original dimensions reap- 
peared at z = 25 mm, gaining power as it propagated. 
The profiles measured in the case of "free- filament" prop- 
agation are reported in Fig. \\j^c) for comparison, in the 
same z range. Note how the effect of the beam stop- 
per is barely detectable after only 20 mm of propagation. 
These results unequivocally show that the observed fila- 
ment can by no means be described in terms of a soliton, 
beam self-guiding effect. 




FIG. 2: Transient stage of the filament formation: (top) ex- 
perimental, P p = 17 P cr ; (bottom) numerical, P = 25 P C r, 
/? (4) ~2.0x 10- 34 cm 5 /W 3 . 

In what follows, we intend to demonstrate how the key 
features of the filament-formation process may be inter- 
preted as the effect of spontaneous evolution of the input 
wave-packet toward a non solit on-like, but still station- 
ary, profile, an hypothesis that has never before been 
considered. In order to make the point clear, we propose 
here the simplest possible model that fully supports such 
a scenario, by adopting the quite severe CW approxima- 
tion, though a priori exclusion of a possible contribution 
of temporal or spatio-temporal (ST) effects was not possi- 
ble. In fact, the results shown in Fig.H are a ls° compati- 
ble with standard, moving-focus |l4| or dynamic spatial- 
replenishment |l5| models (this can be understood by 
noting that if different temporal slices are focused at dif- 
ferent planes, pinholes or stoppers could have different ef- 
fects on different temporal portions of the wave-packet). 
Indeed, close inspection of ST wave-packet profiles by 
means of a high-resolution (20 fs) 3D-mapping technique 
(in the case of input pulses shorter by a factor two) re- 
vealed a dynamics that requires the full 3D model to be 
described precisely [H3- However, such ST effects, the 
importance of which dramatically increases on shorten- 
ing the pulses, are not necessary for filament formation. 
From the experimental viewpoint, this claim is strongly 
supported by the measurements with five times longer 
pulses (we do not present the data here), whose dynam- 
ics appeared to be very similar to the present case, espe- 
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cially in the asymptotic regime. Further support comes 
a posteriori from the rather good matching achieved be- 
tween experiment and CW calculations. As a second, 
highly qualifying, approximation, we have neglected the 
role of the plasma- defocusing effect, thus avoiding the oc- 
currence of any self-guiding regime. In so doing, we wish 
to focus attention on the final stage of the dynamics, 
where plasma-induced effects are supposed to play only 
a minor role, seeking the occurrence of a "final state" 
that behaves as an "attr actor" for the entire dynamics. 
The only non-linear terms included in the model are, 
therefore, those of self-focusing and MPA. Within the 
framework of the paraxial approximation adopted, the 
resulting CW, elliptical, modified non-linear Schrodinger 
equation for the field amplitude A reads: 

dA_J_ fd^A d 2 A \ 
dz 2k \dx 2 dy 2 J 



(i) 



c ■ ' 2 ' ' 

where z is propagation distance, — 2.7 x 10 -16 cm 2 /W 
is the non-linear refractive index, and is the MPA 
coefficient. In our model we took K = 4 in order 
to account for three-photon absorbtion (water band-gap 
W g = 6.5 eV, photon energy 2.4 eV) and, to a first ap- 
proximation, for further absorption due to the photo- 
induced plasma. We verified that similar results are 
also obtained with different K. Owing to the over- 
simplification of our model, which does not allows for 
ah initio evaluation of parameters, for the absorption co- 
efficient we took the value that better fitted the ex- 
perimental results (see the figure caption). The fit led 
to a slight dependence of (3^ on input beam power, 
which may be understood as owing to the different role 
played by plasma absorption. We consider evolution of 
a linearly polarized beam propagating with central fre- 
quency uoq and wave- number k = uoqu/c, with n = 1.334. 
The resulting critical power for CW beam collapse is 
P cr = 3.77A 2 /(87rnn 2 ) = 1.15 MW. The input beam pro- 
file was taken as Gaussian with a 0.12 mm diameter at 
FWHM. Note the slightly larger diameter (and also larger 
powers, see the figure captions) that had to be taken with 
respect to the experiment, in order to compensate for a 
slight overestimation of the losses by our model. Equa- 
tion JU was solved numerically by means of the FFT 
Runge-Kutta split-step method. 

A very surprising result is that our model, designed for 
a qualitative description of the asymptotics, also turned 
out to describe the very transient regime rather well, 
where plasma is eventually supposed to play its major 
role. Both the experimental and numerical results shown 
in Fig. indicate the occurrence of an appreciable flat- 
tening of the top-beam profile in the first 7 mm of prop- 
agation (left column), which occurs in conjunction with 
a slight beam self-focusing. Owing to the flattop pro- 
file, the major effect of this self-focusing becomes that of 



developing a sharp and rather intense ring peak on the 
external border of the "mesa" (see the "volcano" pro- 
file shown in the central column of the figure). Owing 
to the Kerr non-linear response, this ring-like modula- 
tion gives rise to a short, positive, toroidal lens, which 
focuses the outside corona (the "skin") of the beam and 
leads the very intense spike to appear in the center (right 
column in the figure). Owing to its very high intensity, 
this spike leads to a very strong, though only local, MPA. 
The persisting toroidal lens, however, keeps refilling the 
spike, thus establishing a quasi stationary regime (fila- 
ment), which remains as long as most of the beam power 
is coupled to the central spike. 

Figure Efa) presents a comprehensive summary of 
the measured and calculated dynamics, by plotting the 
FWHM beam diameters in the entire (i.e., transient 
plus quasi-stationary) range investigated. Both the free- 
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FIG. 3: (a) HWHM beam radius vs. z: experiment (full cir- 
cles) and simulation (dashed line). Open circles and solid 
line: the same in the case of a 55 u-m pinhole only trans- 
mitting the central spike, (b) Calculated transverse intensity 
profiles for the clipped (left), free (center) and stopped (right) 
filament case. Experiment: P p — 13 P cr - Numerical values: 
P = 15 P cr ; A) ~ 1.25 x 10" 34 cm 5 /W 3 . 

filament and the 55 jam clipped filament are considered. 
Note the rather impressive capability of our model in de- 
scribing the apparent stationary-filament regime, in the 
absence of any defocusing effect. Note also how the model 
predicts quenching of the filament, when the pinhole is in- 
serted. For a more complete description of the dynamics, 
the corresponding calculated transverse intensity profiles 
vs. z are reported in Fig. EJb). Note that here the case 



4 



of the 55 iLtm beam stopper is also depicted (right). 

In order to elaborate the physical interpretation of the 
entire dynamics, in Fig. 0f a) we plot the (measured and 
calculated) peak fluence F p and intensity I pi and the (cal- 
culated only) fractional power losses dP/P vs. z. The 
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effect) and non-linear losses, which continue to absorb 
the generated spike. We interpret the entire dynamics as 
the attempt of the wave-packet to readjust its shape, due 
to the non-linear coupling, according to the requirements 
of maximum localization, minimum losses and stationar- 
ity. We predict the existence of an exact, infinite-power, 
Bessel-like stationary solution of the modified CW non- 
linear Schrodinger equation in the presence of non-linear 
losses. 

The authors acknowledge assistance in numerical cal- 
culations by A. Parola, M. Porras and A. Berzanskis, dis- 
cussions with R. Danielius and A. Piskarskas, technical 
support in the measurements by E. Kucinskas and finan- 
cial support by MIUR (CofinOl/FIRBOl) and EC CEBI- 
OLA (ICA1-CT-2000- 70027) contracts. 



FIG. 4: Calculated peak intensity I p (solid line); measured 
peak fluences F p (full circles); calculated fractional power 
losses (dashed curve), for the same conditions as in Fig. 

results clearly indicate the occurrence of two different 
regimes: for z < 12 mm, both I p (F p ) and dP/P in- 
crease due to the overall effect of self-focusing and of 
MPA. After that, however, while I p (F p ) keeps a high 
and almost constant value, a sharp drop of dP/P occurs. 
Close inspection of the results in Fig. Efb) indicates that 
the beam has undergone an important transformation: 
from a Gaussian to Bessel-like profile, which has pre- 
served appreciable localization (e.g., high I p ) and sta- 
tionarity while minimizing the non-linear losses (MPA). 
In fact, owing to the presence of "cold", slowly decay- 
ing tails, which contain most of beam power, a conical 
wave is provided with a large reservoir that can refuel the 
"hot" central spike and so preserve stationarity, in spite 
of the presence of non-linear losses. If the robustness 
of the Bessel profile under Kerr-induced spatial phase 
modulation is accounted for one could foresee that 
there might indeed exist a Bessel-like, infinite-power, ex- 
act stationary solution of Eq. (Q that behaves as strong 
attractor for the entire non-linear beam transformation. 

In conclusion, by clipping or stopping a light filament 
while it propagates in water we have shown experimen- 
tally that it does not behave as a self-guided wave packet, 
being structurally sustained by a strong energy flux from 
the surrounding beam. The matching obtained between 
the experimental results and those based on a CW model 
that only accounts for self-focusing and non-linear losses 
provides a strong indication that all the temporal effects 
as well as all those related to plasma-induced defocus- 
ing are not essential to the occurrence of the apparent 
guiding effect. In contrast, the quasi-stationary filament 
appears as the dynamical balance between sharp focusing 
of the very external periphery of the beam (skin focusing 
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